7.6.7 Listing (vervolg) 

fCSTBiN IV G LFVEL 1, HO0 3 DDE! 


CCC1 


SUEICUTINF. CD EX ( A, N, NDIK, P, D) 

COO 2 


0 FAL* 8 A (NDIK,NDIM) ,D 

CC 03 


INTEGER P(1) , H 

CCC4 


REAL*8 W D , 5 

C0C5 


REAL W , R , V (256) 

CCC6 


INTEGER £K 

CCC7 


IF (N) 001,001,002 

CCC8 

('01 

WRITE (3,900) 

CC 0 9 

9< 0 

FORM AT (* 1K1 OF N> NDIM OF N>2 56 

CC10 


SIC F 

0011 

0 C 2 

IF (N-NDIH) 003,003,001 

C C 12 

003 

IF (N-256) 004,004,001 

OC 12 

OC 4 

CONTINUE 

C C14 


DC 102 1=1,N 

00 15 


W = C-C 

C C 16 


DC 101 0=1,N 

OC 17 

1C1 

W = S+A (I,J) *A (I,J) 

C C 18 

102 

V (I) = 1.0/SQRT (S)_ 

OC 19 


D= 1 • 0 

CC2C 


DC 111 K= 1 , N 

0C2 1 


R =- 1. C 

CC22 


K1= K-1 

00 23 


DO 105 1=K,N 

CC24 


IF {K— 1) 1031, 1031, 1021 

0025 

102 1 

CONTINUE 

C C 26 


WD=0.0 

0027 


DO 103 J=1,K1 

C C2 £ 

103 

HD=W D+A (I, J)*A<J,K) 

0029 


A (1,K) =A (I, K) - WD 

C C3f 

1031 

CCNTINUE 

0031 


W=A(I, K) 

CC32 


W= AES (W) *V (I) 

0033 


IF (N-R) 105,105,1 04 

l C 34 

104 

F = W 

0035 


P K = I 

C 036 

105 

CCNTINOE 

0 0 37 


P (K) =PK 

CC38 


V (E K ) = V (K ) 

0 0 39 


S=A (PK , K) 

CC4C 


C = 5* D 

0041 


IF (PK-K) 106,108,106 

CC42 

106 

C = - D, 

0042 


DO 107 J =1 , N 

C 04 4 


HD-A (K, J) 

0045 


A (K , J) =A (PK, J) 

CC46 

107 

A (EF, J) = WD 

0047 

1C 8 

KP 1 =K+1 

C C4£ 


IF (F-N) 1081,1101,1101 


It 


D AT L = 


O 
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7.6.7 Listing (vervolg) 

SCblhAH IV G I £VEL 1, MOD 3 


DDET 


C049 

1C 6 1 

CONTINUE 

CC5C 


DC 110 l=Kt1,N 

00 51 


*1D = C . 0 

CC52 


IF IK-1) 1091, 1091, 1002 

0C 5 i 

1C 82 

CONTINUE 

CC54 


CC 109 J=1,K 1 

0055 

109 

fcD = WD+A(K,J)*A (J,I) 

C C 56 

1C 01 

CONTINUE 

0057 

110 

A <E,I) =(A (K,I)-fc'D)/S 

C C56 

1101 

CC NTIKUE 

CC59 

111 

CC NTIK UE 

CC6C 


EETUEK 

0C6 1 


END 


DATE 
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FOR-NUM 38 


70107 

















7.8.7 Listing (vervolg) 

FORTRAN IV G I EVEL 1, M-JiJ 3 DTNV DATE = 70107 


0G0 1 


SU8R0UTI NE D I N V ( A ,N,ND1 R ,P) 

COO? 


REAL*8 A ( NO I M , NO 1 M) 

0002 


INTEGER P(l) 

C004 


REAL*8 WD »R » VI 2 56) 

COO 5 


DO 3 C 6 K I N V = 1 » N 

CC06 


K =N + 1- K I N V 

0007 


K1=K +1 

GCC 8 


T h <K-N) 300,3021 ,3021 

0010 

3U0 

CONTINUE 

C C 1 c 


01 3 C 2 JI N V=K1,N 

0011 


J=K1+N-JINV 

CC 12 


A( J f K 1 ) - V{ J) 

00 13 


WD = C. 0 

GO 14 


00 301 L = K1 ,N 

0015 

301 

WD=A( K,L) *A( L, J) + WO 

GO 16 

302 

V( J)=-WD 

OC 17 

2C21 

CONTINUE 

0018 


R=A(K,K) 

OC 19 


IF ! K-N ) 3022,3041 ,3041 

GO 20 

3C22 

CONTINUE 

002 1 


DO 3C4 JINV=K1,N 

002 2 


J=K 1 + N- JI NV 

CO 23 


A ( Kf J ) - V( J) 

00 24 


WD = C . 0 

00 2 5 


DO 303 L=K1,N 

0026 

303 

WD=A I J , L) *A( L,K) + WD 

0027 

304 

V( J ) =- WD/R 

0028 

2C41 

CONTINUE 

C C 29 


WD = C • C 

CC 30 


IF (K-N) 3042,3051 ,3 051 

003 1 

3C4 2 

CONTINUE 

0 G 2 2 


DO 3C5 L = K1 ,N 

003 3 

3C5 

WD=A( K,L)*A(L,K) + WD 

C 03 4 

3051 

CONTINUE 

00 3 5 

3 C 6 

V( <)=( 1.0-WD) /R 

CO 36 


DO 3C7 JINV=1»N 

CO 37 


J = 1+N-JINV 

C C 38 

307 

A ( J, 1 ) = V( J) 

CO 39 


NM1=N-I 

C04C 


DO 3 1C KINV=1,NM1 

0041 


K=N-KI N V 

CO4 2 


K1=P( K ) 

004 3 


IF (Kl-K) 308,310,308 

0 0 4 4 

3C 9 

DO 3C9 L = 1 , N 

CO 4 5 


WD=A( L,K) 

0C46 


A(L,K)=A(L,K1) 

0 C 4 7 

2C9 

A ( L , K 1 ) = WD 

C 0 4 o 

310 

CONTINUE 
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FGRTRAN IV G LEVEL I, MOD 3 

CC49 RETURN 

CG50 END 


DINV 


DATE = 70107 
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